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Electronic and lattice properties of MgB2 under uniaxial and hydrostatic compression are cal- 
culated. Lattice properties are optimized automatically by using the first-principles molecular 
dynamics (FPMD) method. Features of the electronic band structures under uniaxial and hy- 
drostatic compression are quite different each other. 
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§1. Introduction 

Very recently, a new high T c intermetallic supercon- 
ductor material MgB 2 (T c — 39 K) was discovered by 
J. Akimitsu group.EP MgB2 has a AIB2 crystal sym- 
metry (P6/mmm) and the number of atoms in a unit 
cell is jdhcee. Already, a variety, xif extended experi- 
mentafl&S&B an^heoreticaifflBBBBrJ studies 
have been devoted.lHiP The electronic structures of MgBa 
under high pressiire (hydrostatic)H and with different 
lattice constants^ are calculated by FLAPW, respec- 
tively. 

In this study, we calculate the electronic and lattice 
properties of MgB2 under uniaxial compression by us- 
ing the first-principles molecular dynamics (FPMD) and 
to compare with results under hydrostatic compression. 
It is possible to obtain the optimized lattice properties 
under high pressure conditions using FPMD. We inves- 
tigate the changes of electronic band structures under 
a variety of compression (uniaxial, hydrostatic, varying 
the c/a ratio). Phonon frequencies at the T point under 
uniaxial and hydrostatic compression are calculated in 
this study. 

§2. Method of Calculation 

The present calculation is based on the local .density 
approximation in the density functional theoryOEZP with 
the Wigner interpolation formulaEJ for the exchange- 
correlation. The optimized pseiidopotential (Mg and B) 
by Troullier and Martins (TM)a' is used in order to re- 
duce the number of plane waves. Non-local parts of pseu- 
dopotential are transformed to the Kleinman-Bylander 
separable forms.Ef No ghost bands by this treatment 
are fouiid in our calculation. A partial core correction 
(PCC)EP is considered for Mg pseudopotential. The 
wave function is expanded in plane waves, and the energy 
cutoff is 81 Ry with the maximum number of plane waves 
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being about 2000. The number of sampling k-points is 
95 in the 1/24 of the Brillouin zone (BZ) for all calcu- 
lated systems. In the cell optimization, the number of 
k-points is fixed because a change of the c/a ratio under 
compression is not so large (c/a = 0.868 ~ 1.123). 

The outline of our calculational process with the cell 
optimization is as follows. In this calculation, the system 
has kept to the AIB2 crystal symmetry in the process of 
the cell optimization. Therefore, it does not discuss a 
structural phase transformation under uniaxial and hy- 
drostatic compression in this study. At first, an elec- 
tronic part is quenched to the ground state of a given 
initial lattice parameter, and optimization of electronic 
and unit cell shape parts-is carried out simultaneously 
in the next step. Stresses!^ acting on unit cell surfaces 
are calculated, and the lattice parameters are tuned by 
using them. We can estimate that a numerical error of 
the calculated pressure is less than about 1 GPa. There 
is no effect to discuss qualitatively although this value (1 
GPa) may be slightly large. 

We calculate in three conditions as follows. The first 
is the uniaxial compression. The uniaxial compression 
acts along a direction of c-axis. Values of external com- 
pression along a- and 6-axis are GPa in this case. The 
second is the hydrostatic compression. The third is to 
vary the c/a ratio when the values of a and b are fixed in 
the equilibrium lattice parameters. In the case of varying 
the c/a ratio, all surfaces of the unit cell have non-zero 
internal pressures. 

§3. Results and Discussion 

The optimized lattice properties in this calculation are 
tabulated in Table I. 

Calculated bulk modulus of MgEbJs 145 GPa, which 
agrees with other theoretical result O A heat of forma- 
tion of MgB 2 is 0.61 eV/cell. The calculated equilibrium 
lattice constant of c-axis is shorter by 3 % than that 
of experiment .Ef The number of sampling k-points (= 
95) may be insufficient for discussing quantitatively the 
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Table I. Equilibrium lattice constants[A] and the c/a ratio of 
MgB2. "Pz" indicates the uniaxial compression along c-axis. 
"P" indicates the hydrostatic compression, "c/a" indicates the 
varying c/a. The value of compression "c/a" indicates that along 
c-axis. 
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GPa 
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10 GPa(Pz) 


3.307 


3.053 


1.083 


20 GPa(Pz) 
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0.985 
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20 GPa(P) 


3.246 


2.958 
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30 GPa(P) 


3.179 


2.923 


1.088 


40 GPa(P) 


3.122 


2.894 


1.079 


50 GPa(P) 


3.072 


2.868 


1.071 


22.3 GPa(c/a) 


3.175 


3.047 


1.042 


48.4 GPa(c/a) 


2.910 


3.047 


0.955 


84.6 GPa 
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2.646 


3.047 


0.868 
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3.524 


3.086 


1.142 



electronic and lattice properties because MgB2 is metal- 
lic. In addition, this number is not adjusted in varying 
the unit cell volume. There is no problem for discussing 
qualitatively the electronic band structure changes un- 
der uniaxial and hydrostatic compression in this study. 
To check the convergence, we also calculate in the con- 
dition of the 259 k-points and 100 Ry under GPa. A 
deviation of the equilibrium lattice constants is less than 
about 0.01 A, and slightly improve the error of the lattice 
constant of c-axis from 3 % to 2.6 %. 

The electronic band structures of the optimized lattice 
parameters are calculated. The electronic band struc- 
tures under uniaxial compression (0, 20, 50 GPa) are 
shown in Fig. 1. The electronic band structures under 
hydrostatic compression (20, 50 GPa) are also shown in 
Fig. 1. The electronic band structure (c/a = 0.955, Pz 
= 48.4 GPa) in varying the c/a ratio is also shown in 
Fig. 1. As a result, the electronic band structures at 
20, 50 GPa under hydrostatic compression are similar to 
that at GPa. The small change of the c/a ratio un- 
der hydrostatic compression indicates the isotropic com- 
pressibility. This trend agrees with the other theoretical 
studjlii although the c/a ratio (= 1.123, P = GPa) 
at present study is slightly-different from the other the- 
oretical result (= 1.1468)y ) The widths of band at 20, 
50 GPa is larger than that at GPa (that of 50 GPa 
is largest) as shown in Table II. On the other hand, 
the electronic structures under uniaxial compression are 
quite different from those under hydrostatic compression 
at the Fermi level. The a bands at the T- A line shift 
to lower energy as the pressure increase. The electronic 
structure with varying the c/a ratio (c/a — 0.955, Pz = 
48.4 GPa) is also different from those under hydrostatic 



compression as shown in Fig. 1. This band structure (a 
= an and c = 0.85co, an and c G are the calculated equi- 
librium lattice constant in this study.) is corresponding 
to that of a = an and c = 0.8cotiP (an and Cnl23 ) ) and 
they are similar to each other. 

We calculate density of states (DOS) at GPa, 
50 GPa(hydrostatic), 50 GPa(uniaxial), nearly 50 
GPa(varying c/a, Pz = 48.4 GPa). The density of states 
at Fermi level (N(e/)) at GPa, 50 GPa(hydrostatic), 
50 GPa(uniaxial), nearly 50 GPa(varying c/a) are 1.0, 
0.94, 0.71, 0.77, respectively, where the values of N(e/) 
are normalized as 1.0 for N(e/) at GPa. The number 
of k-points is not enough to discuss quantitatively the 
values of N(e/). The values of N(e/) under compression 
are smaller than that under GPa and the decrease of 
N(ej) under anisotropic compression is larger than that 
under hydrostatic (isotropic) compression. The density 
of states at GPa and 50 GPa(hydrostatic) are similar 
to each other except for the band width. The density of 
states at 50 GPa (uniaxial) and nearly 50 GPa(varying 
c/a) are quite different from those at GPa and 50 
GPa(hydrostatic). 

Table II. The values of the band width from the bottom of the 
valence band to the Fermi level. "Pz" indicates the uniaxial com- 
pression along c-axis. "P" indicates the hydrostatic compression, 
"c/a" indicates the varying c/a. 





Width(eV) 


c/a 


GPa 


12.6 


1.123 


20 GPa(Pz) 


12.9 


1.038 


50 GPa(Pz) 


14.0 


0.908 


20 GPa(P) 


13.3 


1.097 


50 GPa(P) 


14.0 


1.071 


48.4 GPa(c/a) 


13.9 


0.955 



We calculate phonon frequencies / w at the T point un- 
der GPa, 50 GPa(uniaxial) and 50 GPa(hydrostatic) 
compression as shown in Table III. The calculated dis- 
placements of atoms in the unit cell are shown in Fig. ||. 
The number of k-points is 512 (8x8x8) in the whole Bril- 
louin zone (BZ) in order to obtain the accurate phonon 
frequencies because this calculation does not consider 
symmetry due to the displacement of atoms in the unit 
cell. The values of the phonon frequencies at the T 
point. M£,in good agreement with other theoretical re- 
sults&BEJ except for the E2 g mode. The discrepancies 
of the phonon frequencies at the Ei„, A2 U , Bj-jmodes 
at present and the other theoretical studiesDBIiiJ are 
very small. On the other. hand_the values of the F,2 g 
mode (629 [Present], 470,EP 585,0 665§ cm" 1 ) in the 
theoretical calculations has a large difference. The val- 
ues of phonon frequencies under compression are larger 
than that under GPa because the cell volume is shrunk 
by compression. Values of the volume ratio (V/Vn, 
Vo equilibrium volume) are 0.795 (P — 50 GPa) and 
0.872 (Pz = 50 GPa), respectively. The deviations of 
lattice changes are -5.9 % (P = 50 GPa), +2.6 % (Pz = 



Electronic structures of MgB2 under uniaxial compression 



3 



50 GPa) along a-axis and -17.0 % (P = 50 GPa), -10.2 
% (Pz = 50 GPa). The value of the phonon frequency at 
the Fi2 g mode under hydrostatic compression (50 GPa) is 
886 cm -1 , in which the increase of the phonon frequency 
is largest in all calculated phonon frequencies. In partic- 
ular, the change of the phonon frequency at the E 2g mode 
under uniaxial compression (Pz = 50 GPa) is largest in 
the four modes in spite of the elongation of a B-B dis- 
tance on the B layer. The a bands at the T - A line are 
completely occupied by uniaxial compression. This oc- 
cupation means that the number of electrons with regard 
to a bonding in the B layer increases. Increasing a bond- 
ing electrons under uniaxial compression is to strengthen 
the B-B bond, which leads to increase the value of the 
phonon frequency at E2 S mode. On the other hand, the 
lattice change (-5.9 %) under hydrostatic compression 
(P = 50 GPa) directly strengthens the B-B bonding al- 
though the a bands at the T - A line are unoccupied. 



Table III. The values of the phonon frequencies / w (cm -1 ) at the 
r point. "Pz" indicates the uniaxial compression along c-axis. 
"F" indicates the hydrostatic compression. 





/ w (0 GPa) 


/ w (P = 50 GPa) 


/ w (Pz = 50 GPa) 




348 


500 


513 


A 2u 


398 


594 


517 


E 29 


629 


886 


851 


Big 


707 


814 


747 



§4. Summary 

We calculate the electronic and lattice properties of 
MgB2 under uniaxial, hydrostatic compression and vary- 
ing the c/a ratio by using FPMD. The phonon frequen- 
cies at the T point under GPa, 50 GPa(uniaxial) and 50 
GPa(hydrostatic) are calculated. It is possible to explain 
the increasing the phonon frequencies under compres- 
sion by the analysis of the electronic band structures. 
Although the width of the electronic band structure is 
spread by compression, the band structure change is not 
so large under hydrostatic compression. The electronic 
band structure changes (a bands) are quite large under 
uniaxial compression and varying the c/a ratio at the 
Fermi level. We think that it is very important to in- 
vestigate the difference of the electronic band structures 
(tr bands) at the Fermi level under hydrostatic, uniaxial 
and varying the c/a ratio. 

More accurate study (the number of k-points, possi- 
bility of a pressure induced phase transition, etc.) is a 
near future task. 
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MgB2(AlB2) c/a =0.955 MgB2(AlB2) P = 20 GPa MgB2 (A1B2) P = 50 GPa 



Fig. 1. Energy band structures of MgB2 under uniaxial and hydrostatic compression, c/a = 0.955 (varying the c/a ratio). The Fermi 
level is indicated by the horizontal line. 



Electronic structures of MgB2 under uniaxial compression 



5 





Fig. 2. Atomic displacements in the unit cell for phonon frequencies at 
Eiu, A2u, E2g, B13 modes. The arrows indicate the directions of atomic 
displacements in the unit cell. 



